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An association between anisotropic plasma heating and instabilities in the solar wind 
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We present an analysis ol the components of solar wind proton temperature perpendicular and 
parallel to the local magnetic field as a function of proximity to plasma instability thresholds. 
We find that T± p is enhanced near the mirror instability threshold and Tj| p is enhanced near the 
firehose instability threshold. The increase in T± p is consistent with cyclotron-resonant heating, but 
no similar explanation for hot plasma near the firehose limit is known. One possible explanation is 
that the firehose instability acts to convert bulk energy into thermal energy in the expanding solar 
wind, a result with significant implications for magnetized astrophysical plasma in general. 

PACS numbers: 96.60.Vg, 96.50.Tf, 96.50.Ci, 95.30.Qd 



Introduction. — Particle velocity distribution functions 
in the solar corona and solar wind are anisotropic, with 
separate temperatures Xj_ and Tji relative to the mag- 
netic field B [l| . Characterizing the processes that create 
and limit temperature anisotropy is important for under- 
standing heating and dynamical effects in solar physics Q 
and in astrophysical plasmas in general [3| . Comprehen- 
sive in situ measurements of the solar wind by spacecraft 
are a unique way to investigate anisotropic plasmas and 
provide observational constraints for more exotic astro- 
physical objects such as accretion disks [1] and blazars ||. 
The range of anisotropy seen in the solar wind is a result 
of competing phenomena: adiabatic expansion, heating 
through the anisotropic dissipation of waves, instabili- 
ties, and Coulomb collisions. Each of these factors has 
individually been the subject of intensive study; however, 
little work has been done on how the effects interact with 
and regulate each other. Our purpose here is to present 
the first observations of these interactions by looking at 
heating in the presence of instabilities. We begin with a 
review of anisotropy in the solar wind. 

If solar wind protons expanded adiabatically, they 
would conserve the first and second CGL invariants so 
that T± p cx B and T\\ p oc n 2 /B 2 @, where p denotes 
protons and n p is the number density. We would then 
expect R p = T± p /T\\ p oc B 3 /n p . The wind would then 
evolve along a particular trajectory in (/3ii p , i? p )-space, 
where j3» p — n p kBT\\ p / (B 2 /2^q) is the ratio of the paral- 
lel pressure of protons to the magnetic pressure Q • Ob- 
servations have shown that while R p decreases and /3|| p 
increases with distance from the Sun, the slope of the tra- 
jectory is inconsistent with adiabatic expansion. Specifi- 
cally, an additional source of perpendicular heating must 
be introduced 

Evidence for perpendicular heating is common in re- 
gions such as the ionosphere, the solar wind, and the 
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FIG. 1: The distribution function / of the solar wind in 
(J3\\ p , -R P )-space at 1 AU. Lines are theoretical curves of con- 
stant growth rate for the firehose (dotted), mirror (dashed), 
and cyclotron (dot-dashed) instabilities. These instabilities 
bound a stable range of R p as a function of [3« p . 



solar corona, where ions may develop R > 20 One 
explanation is Alfven-cyclotron dissipation, in which ions 
enter cyclotron resonance with compressive Alfven waves 
on spatial scales near their gyro-radius and a re p refer- 
entially energized perpendicular to B 
heating consistent with an Alfven-cyclotron mechanism 
is also directly observed in interplanetary space (l2l |. 

Instabilities prevent expansion and heating from driv- 
ing R p arbitrarily far from unity. Fig. Q] is the probability 
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distribution / of the solar wind in (j3\\ p , i? p )-space as mea- 
sured by the Wind spacecraft at 1 AU. These observa- 
tions are the focus of this letter and are described in more 
detail in the next section. Given the large anisotropics 
in the corona and the effects of adiabatic expansion, one 
would expect R p to vary over orders of magnitude and be 
as small as 10~ 3 @,[l3j]. Instead we find 0.1 < R p < 5 and 
constrained to a narrowing region near isotropy as /3|| p 
increases. This narrowing of / with (3\\ p is attributed to 
instabilities, driven by R pi that generate electromagnetic 
fluctuations, scatter particles in velocity space, and drive 
R p toward unity. For R p < 1 and (3\\ p > 1 the firehose 
instability can limit anisotropy (l4| . while for R p > 1, the 
mirror and cyclotron instabilities are active 0i. One 
way to quantify instabilities is to calculate the rate of 
growth of the unstable modes of the linear Vlasov equa- 
tion. The three curves in Fig. [T]are contours of constant 
growth rate 10 -3 times the proton cyclotron frequency 
for each instability [la ]. The sharp drop in / beyond 
these curves is interpreted as evidence of the instabil- 
ities. Comparing the shape of / with these curves, the 
mirror instability appears to be more important than the 
cyclotron instability at limiting R p > 1, even though the 
amplitude of the cyclotron instability grows more quickly 
for /3|| p < 2. This diminished role for the cyclotron insta- 
bility is supported by the recent discovery of enhanced 
magnetic fluctuations in plasma beyond the mirror and 
firehose threshold curves, but not the cyclotron thresh- 
old [l3j]. This surprising result, which highlights the po- 
tential pitfalls of linear theory, is possibly due to mirror 
fluctuations being more efficient at scattering particles in 
velocity space even if the growth rate is slower [l5j |. 

Having introduced the heating mechanisms and insta- 
bilities associated with temperature anisotropy in the so- 
lar wind, we now explore their interactions. 

Observations and analysis. — Our study is motivated 
by an earlier work that examined the scalar temperature 
T p = (2T± p + T|| p )/3 as a function of R p and (3\\ p and 
produced evidence suggesting plasma near the thresh- 
olds was hotter than expected [13] ■ This would be very 
interesting because the instabilities should not heat the 
plasma, only isotropize it. However, it was difficult to 
tell if the plasma was truly hotter near the instability 
thresholds, or if the result was due to the dependence 
of (3u p on T|| p . To better understand these associations, 
we use a more accurate set of temperature measurements 
and examine the components T± p and Ty p separately. 

This study makes use of observations from two instru- 
ments on the Wind spacecraft: 92-second cadence ion 
velocity spectra from the Solar Wind Experiment Fara- 
day Cup (FC) instruments [18| and 3-second measure- 
ments of B from the Magnetic Field Investigation [l9j |. 
This merged dataset is publicly available and has been 
described in detail elsewhere [20|. In all of our earlier 
work, and to the best of our knowledge in all other de- 
terminations of anisotropy in the solar wind, higher time 
resolution B values were averaged over each ion mea- 
surement. R p was then determined by examining how T p 




FIG. 2: Temperature T p = (2T ±P + Tj| p )/3 as a function of 
Rp and /3u p . Curves indicate theoretical thresholds due to 
instabilities. There is an overall trend of increasing T p with 
/3|| p , which is expected since /3|| p is proportional to Tj| p , but 
there is also a clear association between the thresholds and 
hot plasma. 



depends upon direction relative to the average B. We re- 
alized that the strong magnetic fluctuations generated by 
unstable plasma fl3l ] may create errors in R p and rewrote 
our analysis to use the 3-second values of B. This tech- 
nique produces more accurate and often larger R p when 
there are large fluctuations. While a detailed report on 
this method is in preparation, this letter contains our 
first result. 

About 40% of the 4.1 million Wind ion spectra met 
our criteria for use in this study. We required the un- 
certainties in the derived temperatures to be less than 
10%. We only used periods where Wind was far from the 
Earth's bow shock to avoid magnetospheric contamina- 
tion. Finally, we only used observations where Coulomb 
relaxation was not an important factor. For each mea- 
surement, we calculated the Coulomb collisional age A c 
defined as the number of small-angle Coulomb scatter- 
ings the plasma has exp erienced in the time it took to 
reach the spacecraft (l2j . We then required A c < 0.1. 

We divided the selected observations into a 50 x 50 
grid of logarithmically-spaced bins in the (/3|| p , i? p )-plane. 
Within each bin, we calculated the number of observa- 
tions N and the median value of T p . For Fig. [T]we cal- 
culated / by dividing N by the width in /3|| p and R p 
of each bin. Fig. [5] shows T p for all (/3|| p , i? p )-bins with 
N > 50. Beyond a general tendency for T p to grow with 
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FIG. 3: T± p (a) and T\\ p (b) over the (/3u p , i? p )-plane. The increase in T p is mainly in T± p for plasma beyond the mirror 
threshold and in T\\ p for plasma beyond the firehose threshold. 



/3||p, there are clearly two regions with enhanced T p : one 
along the mirror instability threshold and the other along 
that of the firehose instability. Between these regions, 
even at high /3ii p , we see cooler plasma. For comparison, 
T p w 1.8 x 10 5 K at (/3\\ P ,R P ) = (1, 1) but gets as high 
as T p k 4.6 x 10 5 K, or nearly twice as hot, near the 
instability thresholds. 

Fig. [^conclusively establishes that there is a significant 
enhancement in T p near the thresholds for anisotropy- 
driven instabilities, confirming the earlier suggestion of 
this effect [l7|. There is a striking correlation between 
the region of increased T p above the mirror and firehose 
thresholds and the region of enhanced magnetic fluctua- 
tions found by [l3| . This result further confirms the role 
played by instabilities in limiting R p in the solar wind 
and the idea that the mirror instability is more impor- 
tant than the cyclotron instability. 

We have also calculated T± p and Tj| p over the (P\\ p , R P )- 
plane. The results, which are shown in Fig. [3l are dra- 
matic: the heating near the mirror threshold is almost 
entirely _L to B, with T± p « 6.4 x 10 5 K, while the heating 
near the firehose instability is || to B, with Tji p ~ 7.7 x 10 5 
K. Close inspection of Fig.[3]does show an increase, albeit 
smaller, in the other temperature component. 

Consider the high temperature region near the mir- 
ror threshold. In Fig. Q] we plot the median values of 
T± p and Tu p as functions of R p for all observations with 
3 < (3\\p < 30. Here we can clearly see that in addition 
to an increase in T± p for R p > 1 there is also a slight ni- 



fty to heat the plasma directly, these results suggest the 
following interpretation. First, Alfven waves enter cy- 
clotron resonance with the protons, which increases T± p 
and thus R p . Eventually, Rp is sufficiently large that the 
mirror instability sets in and drives the plasma back to- 
ward isotropy. In doing so, some of the energy deposited 
into T± p by anisotropic dissipation is transferred to Tj| p . 

We now turn to the region near the firehose instabil- 
ity. Previous work suggests that the dominant mecha- 
nism producing R p < 1 is that CGL expansion cools T± p 
more quickly than Tji p , at least until the plasma reaches 
the firehose threshold ■ Fig. 2] shows that T± p is only 
slightly cooler for R p < 1 and 3 < /3ii p < 30, while T« p is 
more than four times higher than typical wind, rising to 2 
MK. This result suggests that a process is increasing Tji p . 
One explanation is that there is an as-yet unknown dis- 
sipation mechanism in the solar wind that can raise Tu p . 
There is little theoretical work on || heating in the solar 
wind but other recent observational studies have reported 
cases of unusually large Tj| p [l2T| , A second possibility is 
that the heating arises directly as a result of the plasma 
being driven into the firehose instability by CGL expan- 
sion. Simulations of expanding wind have suggested that 
the firehose instability reduces I^j^ut also produces a 



high energy tail of particles || to B [21|. Perhaps this high 
energy tail subsequently relaxes and heats the plasma. 

The energy for this speculative || heating could come 
from slowing the solar wind down or by modifying the 
expansion of the magnetic field. The solar wind thermal 
crease in Tu p . Since we do not expect the mirror instabil- energy is typically about 1/100 the kinetic energy. So for 
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Proton Anisotropy R p = T ip / T N 

FIG. 4: Median scalar proton temperature T p (black curve) 
and components T± p (diamonds) and T« p (triangles) as a func- 
tion of R p for all observations with 3 < /3u p < 30. For T± p 
and Tu p uncertainty in the mean is indicated with error bars. 



the largest Tn p about 4% of the bulk energy of the wind 
would have to be converted into thermal energy. 

Conclusions. — We have shown that protons with 
anisotropy beyond the mirror and firehose thresholds are 
3 — 4 times hotter than those in typical solar wind. We 
examined T± p and Tji p and found that both components 
are hotter near each instability, but that most of the ad- 
ditional heating is in T± p for R p > 1 and in T» p for 
Rp < 1. This result is interesting because instabilities 
are not understood to heat plasma themselves. The hot 
plasma is therefore either a signature of the interaction 
of instabilities with secondary processes or evidence that 
instabilities do more than merely redistribute thermal 
energy. 

For R p > 1 we suggest that the high temperatures are 
due to a combination of an ion-cyclotron resonant heating 
process (such as dissipation of Alfven waves) increasing 
Tip and redistribution of thermal energy from T± p into 
2]|p by the mirror instability. A plasma undergoing slow 



_L heating may achieve a state where the rate of injection 
of energy into T± p is balanced by the redistribution of 
energy from T± p into T\\ p by the mirror instability. Mea- 
surements of anisotropics alone therefore underestimate 
the level of heating from dissipation. 

For R p < 1 the situation is less clear. Either there is 
an as-yet unidentified parallel heating mechanism at work 
in the solar wind, or the combination of CGL expansion 
and the firehose instability are capable of converting a 
small fraction of the bulk kinetic energy of the wind into 
thermal energy. 

These results have implications beyond the solar wind. 
In any expanding plasma with /3 > 1, the CGL-firchosc 
association may be an effective mechanism for paral- 
lel heating and the generation of magnetic fluctuations. 
This result is directly applicable to heating and magnetic 
field generation in clusters of galaxies [22|]. Our results 
also suggest a heating mechanism for contracting astro- 
physical plasmas. Consider, for example, accretion onto 
a compact object such as a black hole or a neutron star. 
There is a great deal of kinetic energy gained by the ac- 
created matter as it falls into the gravitational potential 
well of the compact object. However, it is not understood 
how this kinetic energy is converted into thermal energy 
in the dynamical time of the infall. We have shown how 
expansion of the solar wind drives the plasma into the 
firehose instability and ultimately heats it. In the case of 
accretion, contraction of the plasma would instead drive 
R p > 1, where the mirror instability might play an anal- 
ogous role to the firehose instability in the solar wind, 
heating the accreting material while slowing the flow or 
shearing the magnetic field. 
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